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Abstract 

A  dynamic  electrochemical  simulation  model  of  a  grid  independent  proton  exchange  membrane  (PEM)  fuel  cell  power  plant  is  presented. 
The  model  includes  the  methanol  reformer,  the  PEM  stack,  and  the  power  conditioning  unit.  The  model  is  then  used  to  predict  the  output 
voltage  and  study  the  transient  response  of  a  PEM  power  plant  when  subjected  to  rapid  changes  in  a  residential  load  connected  to  it.  The 
results  show  the  fast  response  capabilities  of  the  PEM  power  plant  in  following  changes  in  the  load. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  ever  increasing  demand  for  electrical  energy  and  the 
intense  competition  between  electric  companies  in  the  new 
electric  utility  market  has  intensified  research  in  alternative 
sources  of  electrical  energy  that  are  reliable  and  cost  effective. 
The  fuel  cell,  as  a  renewable  energy  source,  is  considered  one 
of  the  most  promising  sources  of  electric  power.  Fuel  cells  are 
not  only  characterized  by  higher  efficiency  than  conventional 
power  plants,  but  they  are  also  environmentally  clean,  have 
extremely  low  emission  of  oxides  of  nitrogen  and  sulfur  and 
have  very  low  noise. 

The  fuel  cell,  in  its  elemental  form,  consumes  hydro¬ 
gen  as  a  fuel  and  produces  dc  power.  The  commonly 
available  fuel  cells  can  be  classified  according  to  tempera¬ 
ture,  into  low-temperature,  medium- temperature,  and  high- 
temperature  fuel  cells  [1].  The  low-temperature  fuel  cells  in¬ 
clude  the  alkaline  fuel  cell  (AFC),  and  solid  polymer  fuel  cell 
(SPFC)  or  proton  exchange  membrane  fuel  cell  (PEMFC). 
The  medium-temperature  class  has  the  phosphoric  acid  fuel 
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cell  (PAFC).  The  high- temperature  class  has  the  molten- 
carbonate  fuel  cell  (MCFC)  and  solid  oxide  fuel  cell  (SOFC). 

The  main  components  of  the  fuel  cell  system  include 
the  fuel  processing  unit  or  the  reformer,  the  fuel  cell  stack 
and  the  power  conditioning  unit.  In  general,  hydrogen  fuel 
is  produced  by  processing  some  hydrocarbon  fuel,  such  as 
propane,  natural  gas,  or  methanol  in  the  reformer.  During 
fuel  processing,  carbon  monoxide  is  produced.  Reduction 
of  carbon  monoxide  to  acceptable  levels  is  achieved  by 
the  water-gas  shift  reaction.  The  low-temperature  and  the 
medium-temperature  fuel  cells  have  external  reformers  for 
producing  hydrogen.  The  high-temperature  fuel  cells  pro¬ 
cess  the  hydrocarbon  fuel  internally  due  to  the  high  working 
temperature.  The  output  from  a  fuel  cell  is  dc  power.  When 
a  fuel  cell  power  plant  (FCPP)  provides  power  to  a  residen¬ 
tial  load,  or  to  the  electrical  grid,  a  power  conditioning  unit 
is  needed.  The  power  conditioning  unit  is  simply  a  dc/dc 
converter  used  to  raise  the  C  output  voltage,  which  is  gen¬ 
erally  the  dc  bus  voltage,  followed  by  a  single-phase  or  a 
three-phase  dc/ac  inverter.  In  grid-parallel  operation  mode,  a 
transformer  is  needed  because  of  the  voltage  difference  be¬ 
tween  the  FCPP  output  voltage  and  the  grid  voltage.  Due  to 
the  low  working  temperature  (80-100  °C)  and  fast  start  up, 
PEMFCs  are  the  best  candidate  for  residential  and  vehicular 
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applications.  The  reasons  behind  selecting  the  PEM  FCPP  in 
this  paper  are:  (a)  fuel  cell  user  groups  publications  indicated 
that  for  small  residential  applications,  where  major  heating 
appliances  are  natural  gas  operated,  a  5  kw  unit  is  adequate, 
(b)  A  5  kw  commercial  PEM  FCPP  is  currently  operational 
in  the  authors’  lab. 

Many  models  have  been  proposed  to  simulate  the  fuel  cell 
in  the  literature.  The  basis  of  a  model  can  be  fluid  dynamics, 
electrochemical  reaction,  heat  transfer  and  thermal  [2-10]. 
To  study  the  transient  response  and  load  following  ability  of 
a  PEM  FCPP,  this  paper  introduces  an  electrochemical  model 
for  a  5  kW  fuel  cell.  The  scheme  incorporates  the  model  for 
an  external  reformer  to  generate  hydrogen  from  methanol. 

The  paper  is  organized  as  follows:  Section  II  introduces 
a  model  for  the  fuel  cell  system.  Load  following  tests  of  the 
model  is  presented  in  Section  III.  Section  IV  presents  the 
conclusions. 


2.  Fuel  cell  system  model 


Fig.  1.  The  PEM  fuel  cell  model. 


where  is  N0 :  number  of  series  fuel  cells  in  the  stack,  I  stack 
current  (A),  F  Faraday’s  constant  (Ckmol-1),  kr  modeling 
constant  (kmol  (s  A)-1). 

Using  equation  (4),  Eq.  (3)  can  be  rewritten  in  the  s  domain 
as: 


2.7.  Fuel  cell  model 


In  [6],  Padulles  et  al.  introduced  a  model  for  the  SOFC. 
This  model  can  be  modified  to  simulate  the  PEM  FCPP  as 
follows: 

The  proportional  relationship  of  the  flow  of  gas  through  a 
valve  with  its  partial  pressure  can  be  stated  as  [6] : 


<7h2 

Ph2 


y/M* 


and 

<7h2q 

Ph2o 


(1) 

(2) 


where  qu2i$  molar  flow  of  hydrogen  (kmols-1),  qu2o  mo¬ 
lar  flow  of  water  (kmols-1),  pn2  hydrogen  partial  pres¬ 
sure  (atm),  pu2o  water  partial  pressure  (atm),  ku2  hydro¬ 
gen  valve  molar  constant  (kmol  (atm  s)-1),  ku2o  water 
valve  molar  constant  (kmol  (atm  s)-1),  kan  anode  valve  con¬ 
stant  (v/KmofKg  (atms)-1),  Mh2  molar  mass  of  hydrogen 
(kg  kmol-1),  Mh2o  molar  mass  of  water  (kg  kmol-1). 

For  hydrogen,  the  derivative  of  the  partial  pressure  can  be 
calculated  using  the  perfect  gas  equation  as  follows  [6] : 


=  (gg  -  2 kr  I)  (5) 

1  +  Tu2S  2 

where: 


Using  Eq.  (5),  the  equations  for  the  partial  pressures  of 
water  pn2o ,  and  oxygen  po2  can  also  be  derived. 

In  [7]  and  [8],  the  authors  introduce  a  model  that  describes 
the  polarization  curves  for  the  PEM  fuel  cell  where  the  fuel 
cell  voltage  is  the  sum  of  three  terms,  the  Nernst  instantaneous 
voltage  E  in  terms  of  gas  molarities,  activation  over  voltage 
77act,  and  ohmic  over  voltage  r/ohmic-  In  mathematical  form, 
polarization  curves  can  be  expressed  by  the  equation: 


Ucell  —  E  -\-  Tjact  T  Johnnie  (2) 

where  rj act  is  a  function  of  the  oxygen  concentration  Co2 
and  stack  current  I  (A),  and  q 0hmic  is  a  function  of  the  stack 
current  and  the  stack  internal  resistance  Rmt  (£2).  Assuming 
constant  temperature  and  oxygen  concentration,  Eq.  (7)  can 
be  rewritten  as  [8]: 

Vceii  =  E  —  Bln(CI)  —  RintI  (8) 


d  RT 

=  (3) 

where  R  is  the  universal  gas  constant  ((1  atm)  (kmol  K)-1),  T 
absolute  temperature  (K),  Van  volume  of  the  anode  (1),  q1^ 

hydrogen  input  flow  (kmols-1),  q hydrogen  output  flow 

(kmol  s-1),  qfni  hydrogen  flow  that  reacts  (kmol  s-1). 

The  relationship  between  the  hydrogen  flow  and  the  stack 
current  can  be  written  as: 


where  constants,  B  =  0.04777  V  and  C  =  0.0136  A  1  [7].  The 
Nernst  voltage  in  terms  of  gas  molarities  can  be  written  as 
[6]: 


PH2P02 

Pu2o 

where  E0  is  the  open  cell  voltage  (V)  and  R  is  the  universal 
gas  constant  (J  (kmol  K)-1 . 

Using  Eqs.  (5),  (8)  and  (9)  the  fuel  cell  model  can  be  drawn 
as  in  Fig.  1  [6]. 
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Fig.  2.  The  reformer  and  the  reformer  controller  model. 


2.2.  Reformer  model 


In  [9]  the  author  introduced  a  simple  model  of  a  reformer 
that  generates  hydrogen  through  reforming  methanol.  The 
model  is  a  second  order  transfer  function.  The  mathematical 
form  of  the  model  can  be  written  as  follows: 


qn2  _  CV 

^methanol  UT2S2  +  (T\  +  T 2)s  +  1 


(10) 


where  ^methanol  is  methanol  flow  rate  (kmol  s-1),  CV  conver¬ 
sion  factor  (kmol  of  hydrogen  kmol-1  of  methanol),  ti,r2 
time  constants  (s). 

To  control  hydrogen  flow  according  to  the  output  power 
from  the  fuel  cell,  a  feedback  from  the  stack  current  is  con¬ 
sidered.  The  relationship  between  the  required  hydrogen  and 
the  increase  in  the  stack  current  can  be  written  as: 


req  N0 1 

qu 2  ~~  2FU 


(11) 


where  q ^  is  amount  of  hydrogen  flow  required  to  meet  the 

load  change  (kmol  s-1),  U  Utilization  rate. 

The  amount  of  hydrogen  required  to  meet  the  load  change 
can  be  used  to  control  the  methanol  flow  rate  [9] : 


^methanol  -  (*3  +  ^  «H2)  (12) 

where  k 3  is  the  PI  gain  and  T3  is  the  time  constant  of  the 
PI  controller.  A  proportional  integral  (PI)  controller  is  used 
to  control  the  flow  rate  of  methanol  in  the  reformer  [9] .  The 
oxygen  flow  is  considered  using  the  hydrogen,  oxygen  flow 
ratio  rh-0.  The  transfer  function  of  the  controller  with  the 
current  feedback  is  given  in  Eq.  (12).  The  reformer  and  the 
reformer  controller  are  illustrated  in  Fig.  2. 


2.3.  The  power  conditioning  unit  model 


x 


Fig.  3.  Fuel  cell,  inverter  and  load  connection  diagram. 


response  accuracy.  On  the  other  hand,  it  complicates  the  sys¬ 
tem  model.  All  commercial  FCPPs  have  to  conform  to  IEEE 
standard  #P- 1547,  which  guarantees  the  ripples  in  the  FCPP’s 
output  voltage  to  be  within  commercially  acceptable  limits. 
A  simple  model  of  the  inverter  is  given  in  [10],  where  output 
voltage  and  output  power  are  controlled  using  the  inverter 
modulation  index  and  the  phase  angle  8  of  the  ac  voltage, 
Vac.  Considering  the  fuel  cell  as  a  source,  the  inverter  and 
load  connection  is  shown  in  Fig.  3.  The  output  voltage  and 
the  output  power  as  a  function  of  the  modulation  index  and 
the  phase  angle  can  be  written  as: 


Vac  —  IP  Vcell  2-8 

(13) 

mV cell  14  . 

Pac  =  sin(S) 

A 

(14) 

Pi 

’  (Vscos  (6»)) 

(15) 

I  =  mlt  cos (6  +  5) 

(16) 

where:  Vac  is  ac  output  voltage  of  the  inverter  (V),  m  inverter 
modulation  index,  8:  phase  angle  of  the  ac  voltage  mVceii 
(rad),  Pac  ac  output  power  from  the  inverter  (W),  Vs  load 
terminal  voltage  (V),  X  reactance  of  the  line  connecting  the 
fuel  cell  and  the  load  (£2),  7t  load  current  (A),  6  load  phase 
angle  (rad),  Pj  load  power  (W). 

As  mentioned  before,  the  modulation  index  is  used  to  con¬ 
trol  the  magnitude  of  the  ac  output  voltage,  which  in  turn  con¬ 
trols  the  amount  of  reactive  power  flow  to  the  load.  The  phase 
angle  is  used  to  control  the  active  power  flow  from  the  in¬ 
verter  to  the  load.  The  voltage  at  load  terminals  is  considered 
constant.  PI  controllers  are  used  to  control  the  modulation 
index  and  the  phase  angle  feedback  signals  to  the  inverter. 
The  transfer  function  of  the  modulation  index  and  the  phase 
shift  are  given  in  Eqs.  (15)  and  (16). 


m  = 


+  kas 

- - ~(Vr 


(15) 


The  power  conditioning  unit  is  used  to  convert  dc  output 
voltage  to  ac.  As  mentioned  before,  the  power  conditioning 
unit  includes  a  dc/dc  converter  to  raise  dc  output  voltage  to 
dc  bus  voltage,  followed  by  a  dc/ac  inverter  to  convert  dc  bus 
voltage  to  ac.  In  this  paper,  only  a  simple  model  of  a  dc/ac 
inverter  is  considered  for  the  following  reasons:  the  dynamic 
time  constant  of  inverters  is  of  the  order  of  microseconds  or 
at  the  most  milliseconds.  The  time  constants  for  the  reformer 
and  stack  are  of  the  order  of  seconds.  Hence,  including  the 
inverter  time  constant  will  have  negligible  effect  on  the  time 


h  +hs/D 

r 

S 


Pac  +  Pi) 


(16) 


where  £5,  k&  £7,  and  k%  are  constants,  Vr  is  the  reference  volt¬ 
age  signal,  Pv  is  the  reference  power  signal,  and  P\  represents 
the  load  power.  The  block  diagram  of  the  inverter  with  the  PI 
controllers  is  illustrated  in  Fig.  4. 

The  current  feedback  signal  in  Fig.  1  and  Fig.  2,  is  cal¬ 
culated  from  the  load  information  and  the  ac  output  voltage. 
Due  to  the  assumed  grid-independent  operation  the  trans¬ 
former  model  is  not  included. 
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Fig.  4.  The  dc/ac  inverter  model. 


3.  Tests  and  Results 

The  models  in  Figs.  1,  2,  and  4  are  connected  in  cascade 
to  form  a  complete  model  for  the  PEM  fuel  cell  system. 

The  PEM  FCPP  we  currently  have  in  the  lab,  as  well  as 
commercial  FCPPs,  need  at  least  an  hour  to  start  from  cold 
to  build  up  the  reformer  and  stack  temperatures.  A  mini¬ 
mum  load  value  (critical  load)  has  to  be  maintained  all  the 
time  during  FCPP  operation.  Decreasing  the  load  below  the 
critical  load  value  will  put  the  FCPP  in  a  sleep  or  dormant 
mode,  where  a  small  amount  of  fuel  is  used  to  maintain  the 
reformer  and  stack  temperatures  at  operating  levels.  In  the 
following  test  cases,  the  FCPP  is  assumed  to  be  in  the  ac¬ 
tive  mode  and  the  initial  active  and  reactive  output  powers 
are  0.8  and  0.6  p.u,  respectively.  The  commercially  available 
PEM  FCPPs  come  equipped  with  storage  batteries  connected 
in  parallel  with  the  dc  bus.  These  batteries  serve  as  a  short- 
period  auxiliary  power  source  to  meet  load  demand  that  can¬ 
not  be  met  by  the  FCPP,  particularly  during  transient  periods. 
Choosing  the  control  system  parameters  affects  the  system 
performance.  For  example,  increasing  the  gain  of  the  PI  con¬ 
troller  that  controls  the  modulation  index  makes  the  voltage 
less  sensitive  to  load  variations.  Testing  the  proposed  model 
indicates  that  the  reformer  and  the  reformer  controller  pa¬ 
rameters  have  a  large  effect  on  the  FCPP  time  response.  The 
reformer  parameters  affect  the  reformer’s  damping  factor. 
Choosing  a  damping  factor  of  1.0  for  the  uncontrolled  re¬ 
former  (open  loop)  will  ensure  that  no  unrealistic  overshoot 
occurs  [9].  This  will  yield  x\  =  T2.  For  maintaining  stabil¬ 
ity  and  reasonable  fuel  processor  dynamics  of  the  controlled 
reformer,  the  controller  parameters  are  chosen  to  ensure  a 
damping  factor  of  0.707  [9].  With  this  assumption,  the  PI 
controller  time  constant  T3  is  equal  to  x\  and  £3  is  equal  to 
(2 CV)-1.  The  model  parameters  are  given  in  Table  1. 

3.1.  Casel 

The  FCPP  model  is  tested  with  step  changes  in  the  load  as 
shown  in  Fig.  5.  These  abrupt  changes  in  active  power  are  for 
testing  the  dynamic  response  of  the  system,  and  do  not  nec¬ 
essarily  represent  changes  in  a  residential  load.  The  change 
of  the  current,  the  voltage,  the  output  power,  and  the  flow  rate 
of  hydrogen  and  methanol  are  illustrated  in  Figs.  6-10.  The 


Table  1 

Model  Parameters 

Stack  Temperature,  (K) 

343 

Faraday’s  constant,  F  (Ckmol-1) 

96484600 

Universal  gas  constant,  R  (Jkmol-1  K) 

8314.47 

No  load  voltage,  E0,  (V) 

0.6 

Number  of  cells,  Na 

88 

Kr  constant  =  N0/4F  (kmol  s-1  A) 

0.996  x  10“6 

Utilization  factor,  U 

0.8 

Hydrogen  valve  constant,  Khi  (kmol  s_  1  atm) 

4.22  X  10"5 

Water  valve  constant,  &h2o  (kmol  s_ 1  atm) 

7.716  x  10“6 

Oxygen  valve  constant,  kcb  (kmol  s-1  atm) 

2.11  x  10"5 

Hydrogen  time  constant,  tho  (s) 

3.37 

Water  time  constant,  th2o  (s) 

18.418 

Oxygen  time  constant,  to  (s) 

6.74 

Reformer  time  constant,  t\  (s) 

2 

Reformer  time  constant,  T2  (s) 

2 

Conversion  Factor,  CV 

2 

Activation  voltage  constant,  B  (A-1) 

0.04777 

Activation  voltage  constant,  C  (V) 

0.0136 

Stack  internal  resistance,  Rint  (£2) 

0.00303 

Line  reactance,  X  (£2) 

0.05 

PI  gain  constants  k$,  ke 

10 

Voltage  reference  signal,  Vr  (p.u) 

1.0 

Methane  reference  signal,  <2methref  (kmols-1) 

0.000015 

Hydrogen-oxygen  flow  ratio,  rh-0 

1.168 
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Fig.  5.  Load  step  changes. 
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Fig.  6.  Fuel  cell  dc  output  current. 
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Fig.  7.  Fuel  cell  dc  output  voltage. 
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Fig.  10.  Methanol  flow  rate. 


Time  [s]  Fig.  11.  Residential  load  profile  for  single  circuit. 


Fig.  8.  ac  output  power. 
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Fig.  12.  ac  output  power. 


Fig.  9.  Hydrogen  flow  rate. 

time  interval  0  to  17  s  is  not  shown  in  the  plots  because  it 
represents  the  dormant  state  of  the  FCPP.  From  these  figures, 
it  is  evident  that  the  increase  in  the  load  increases  the  feed¬ 
back  current,  which  in  turn  decreases  the  output  voltage  of  the 
fuel  cell.  The  increase  in  the  current  increases  the  methanol 
flow  rate  and  the  hydrogen  flow  rate  to  increase  power  flow 
from  the  cell  to  the  load.  As  seen  in  the  ac  output  power 
curve,  the  output  power  has  a  time  delay  in  following  the 
load  power.  This  is  due  to  the  reformer  and  the  fuel  cell  time 


constants. Case  2In  this  test  case,  the  model  is  tested  using  an 
actual  residential  load  profile.  A  load  profile  for  a  2500  sq. 
ft.  house  with  all  electric  appliances  occupied  by  two  adults 
and  four  children  is  used.  Due  to  the  limited  output  power 
of  the  fuel  cell  (5  kW),  only  one  120  V  circuit  of  the  house 
is  considered.  The  active  power  load  profile  for  a  4  h  period 
with  a  15  s  sampling  interval  is  shown  in  Fig.  11.  The  power 
and  voltage  bases  are  5  kW  and  120  V,  respectively. 

The  purpose  of  this  test  is  to  check  the  load-following 
characteristics  of  the  FCPP.  The  output  power,  current  and 
voltage  curves  of  the  fuel  cell  are  shown  in  Figs.  12-14. 
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Fig.  13.  Fuel  cell  dc  output  current. 
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Fig.  14.  Fuel  Cell  dc  output  voltage. 


Comparison  of  Figs.  11  and  12  indicates  that  the  FCPP 
demonstrates  rapid  response  to  residential  load  changes  and 
exhibit  good  load-following  capabilities. 


4.  Conclusions 

This  paper  introduces  a  dynamic  model  for  a  5  kW  PEM 
fuel  cell  system.  The  proposed  dynamic  model  includes  the 
fuel  cell  model,  the  gas  reformer  model,  and  the  power  con¬ 


ditioning  unit.  To  study  the  time  response  of  the  fuel  cell, 
the  proposed  model  has  been  tested  with  step  increase  and 
decrease  of  the  electric  load.  The  obtained  results  show  a  fast 
response  of  the  fuel  cell  to  load  changes.  To  check  the  load¬ 
following  characteristics,  the  proposed  model  has  been  tested 
using  a  typical  residential  load  for  a  period  of  4  h.  The  model 
showes  a  high  degree  of  conformity. 
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